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Steinhausen’s  hypothesis,  that  a  simple  torn  ion -pendu  1  nm  mode]  adequately  describes  the  dynamics  of 
semicircular-canal  action,  has  received  its  share  of  criticism.  This  is  primarily  because  of  the  fact 
that  the  time  courses  of  afferent  discharges  and  vestibular  reactions  to  angul ar-accel erat i on  stimuli 
only  approximately  relate  to  the  time  course  obtained  as  a  solution  to  the  torsion-pendulum  model  for  the 
same  angular  accelerations.  Thus,  there  is  a  need  to  determine  a  generalized  transfer  function  whi'‘h 
delineates  those  components  that  provide  meaning  tw  steinhausen’s  hypothesis  from  those  which  describe 
other  phenomena  such  as  "adaptation". 


t 


INTRODUCTION 

Rotational  head  movements  stimulate  the  sensory  end  organs  or  the  semicircular  canals  of  the  inner  ear, 
causing  compensatory  eye  movements  (nystagmus ) ,  appropriate  to  the  plane  of  rotation,  via  a  vestibulo- 
ocular  reflex  arc.  However,  the  sensory  cpithelia  -r  the  semicircular  canals,  the  cristae  ampul  lares, 
are  relatively  inaccessible  for  direct  functional  .study  during  normal  physiological  head  movements. 

In  the  first  place,  they  lie  deep  within  the  skull,  enclosed  by  Sf*on«y  and  petrous  hone  and  the  osseous 
labyrinth.  Secondly,  any  technique  which  exposes  the  hair  cells  or  their  nerv**  supply  for  neurobiological 
experimentation  necessarily  disrupts  the  vascular  supply  to  the  neuroepithel ium  and  interferes  with  the 
ionic  balance  between  the  perilymph  and  endolymph.  Thirdly,  even  when  the  membranous  ampullae  are 
exposed,  it  is  difficult  to  perceive  the  extent  to  which  the  normally  invisible  cupula  covers  the 
crista  and  fills  the  ampulla  of  the  duct.  (It  is  the  cupula  which  couples  the  mechanical  energy  of  head 
motion  to  the  receptor  hair  cells  of  the  crista  for  transduction  into  neural  impulses  at.  the  afferent 
terminals.)  Thus,  it  is  not  surprising  that,  for  over  a  century,  the  dynamic  response  characteristics 
of  the  semicircular  canals  have  not  been  completely  elucidated  despite  continuous  investigation  by 
many  vestibular  physiologists. 

Historically,  until  the  time  of  Fluorens  the  semicircular  canals  were  associated  with  the 

perception  of  sound.  By  sectioning  the  membranous  semicircular  canals  in  the  pigeon  and  the  rabbit, 

Fluorens  established  that  eye,  head,  and  body  movements  always  occur  in  the  same  plane  as  that  of  the 
injured  canal.  Holtz  (l6),  also  investigating  the  pigeon,  first  associated  these  canals  with  bodily 
equilibrium.  Breuer  (5),  Crum  Brown  (9),  and  Mach  (.’6)  suggested  that  the  crist.ae  ampul  lares 
respond  to  rotatory  accelerations  of  the  head.  Moreover,  they  also  hypothesized  that  the  motion  of 
the  endolymph  within  the  canal  was  responsible  for  eliciting  the  neural  activity  in  the  nerve  endings 
(hydrodynamic  theory).  (Mach,  however,  later  appears  to  have  favoured  a  hydrostatic  theory,  i.e., 
that  pressure  is  the  adequate  stimulus — see  refs.  7°,  Pfi.  )  Fwald  (ll),  using  an  "hydraulic  hammer", 
produced  artificial  endolymphatic  movements  in  the  canals  of  pigeons  and  demonstrated  that  the 
direction  of  head  nystagmus  was  dependent  upon  the  direction  of  endolymph  flow.  Fvald’s  "law"  was 
placed  on  a  firmer  footing  when  Lovenstein  and  Hand  (95)  established  that  the  action-potential 
(impulse)  frequency  in  fibers  Innervating  the  crist.ae  varied  according  to  the  direction  of  angular 
movement  of  the  head. 

In  order  to  describe  the  dynamics  of  the  semiei rcu| ar-canal  system,  fite in hausen  ('6,  ^7),  and, 
subsequently,  F.gmond ,  Oroen  and  .lonkees  (10),  proposed  the  so-called  toralon-pendul  um  mode].  In 
this  model,  the  Instantaneous  angular  deflection  of  the  cupula,  E(V),  mused  by  an  instantaneous 
angular  acceleration  of  the  head,  a(t),  is  governed  by  the  linear  differential  equation 

3  d  |(-t)  t  H  Aiiil  *  AC(t  )  .  fOa(t).  (1) 

dt  dt 

where  0  is  the  effective  moment  of  Inertia  of  the  endolymph  in  the  semicircular  canal;  n  is  the 
viscous  damping  momeni  per  unit  angular  velocity  of  endolymph  relative  to  the  wall  T  the  membranous 
canal;  A  is  the  el ast Ic-reator ing  moment  per  unit  angular  displacement  of  the  end  -lymph  relative  to 
the  membranous  ampul lary  walls;  and  t  is  a  constant  of  proportionality  relating  ratio  volumetric 
displacements  (between  the  cupula  and  the  canal  endolymph).  The  transfer  function  of  Fq ,  (i)  can 
be  written  as : 
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where  A f ;  1  and  “ir  i  «»*•*  IJtpln*.-  *  •*»»/!  ;  ?  rn,s  f  .  t '  *  and  fit),  respectively;  ■’  -  t  ij  t  and  tj  -  n/fl 
and  t  .  -  k>‘|l  hih*  the  ;=  :  fn*  in<i  sh  rt  t  iiv  constants  .f  the  h«*avi  ly-damped  cupul a/end1 -lymph 

sy.it.em  (  |0,  .’Oj, 

•ne  way  t..  t  »•_**  Uie  adequacy  this  m<-del  is  t.>  rec  .rd  the  neural  activity  ,f  the  primary  afferent 
filers  innervut  itif*  the  .,«-micl rculur  canals  a.;  they  respond  t ...  appropriate  stimulation.  This  has  been 
dune  |*<>r  a  variety  •>f  species  >f  animals  I,  12,  17,  IP,  .’3,  2m,  11,  3l< ,  irjt  1ft)  and  the  model  has  been 

found  to  tie  deficient..  The  present  set  <>f  experiments  was  designed  t-  study  the  me chano neural  response 
characteristics  of  primary  afferent,  semicircular-canal  units  in  the  pigeon,  in  order  to  determine 
whether  r  not.  a  general  i/.isi  transfer  function  could  he  obtained  which  would  also  describe  similar 
neurodynamicf.  in  other  species.  (In  testing,  to  this  model,  the  assumption  is  made  that  the  afferent 
activity  is  proportional  to  the  cupular  di splacement . ) 


METHODS 

White  King  pigeons  ( Col umba  ]  ivia )  were  surgically  prepared  for  microelectrode  recording  from 
peripheral  units  which  innervated  the  semicircular  canals  (primarily,  the  anterior  semicircular  canal). 
The  animal  (with  head  immobi  i  ized )  was  oriented  on  board  a  rotatory  device  so  that  the  center  of  its 
head  plane  was  coincident  with  the  plane  of  rotation  (according  to  the  type  of  canal  being  recorded 
from ) . 

The  main  rotatory  sequence  consisted  of  a  scries  of  sinusoidal  angular  accelerations, 

a(t )  *  a  sin  fit  ft ,  (3) 

m 

which  were  delivered  to  anesthetized  preparations  at  frequencies,  f,  from  0.01  to  10  Hz  with  peak 
angular  accelerations,  am,  of  2.0,  !*.0,  8.0,  12.0,  and  20,0  degrees/s^  (t  3  real  time  in  Eq.  C  B ) )  • 

The  single  unit  neural  activity  was  amplified,  displayed  on  an  oscilloscope  and  recorded  on  magnetic 
tape  according  to  conventional  techniques.  (One  channel  of  the  tape  recorder  was  used  for  voice 
commentary;  another  for  the  stimulus  reference  signal.) 

The  stimulus  reference  signal  was  used  to  trigger  a  physiological  signal  analyzer  (Nicolet 
Instruments  Inc.)  to  count  and  store  the  taped,  entrained  action  potentials  in  preselected  appropriate 
time  periods  (bins)  for  up  to  UO96  sequential  bins  (see  Fig.  1  for  typical  binned  response). 

Fourier  techniques  were  used  on  the  binned  neural  data  to  obtain  the  magnitude  of  the  peak  amplitude  of 
the  response  and  the  temporal  (phase)  relationship  between  the  angular  acceleratory  stimulus  and  the 
fundamental  component  of  the  neural  response. 

Amplitude-  and  phase-values  were  usea  in  a  curve-fitting  program  on  a  PDP-11  /b0  minicomputer 
(Digital  Equipment  Corp. )  to  provide  a  distinct  mathematical  expression  for  the  best-fitting  transfer 
function  for  a  linear  system.  (The  steady-state  frequency  response  Tor  a  linear  system  to  a  sinusoidal 
input  angular  acceleration  may  be  found  from  gain, 

!r,(r)|  »  ((He  off))'"'  ♦  dm  i;(f))'’l’>,  (■*) 

and  phase , 

0(f)  ,  tan"1  (Im  <>(f)/Re  0(f)),  (5) 

spectra  ( Hode  plots — see  ref.  30),  where  He  0(f)  and  !m  0(f)  are  the  real  and  imaginary  parts  of  G(f), 
respectively.  ) 

The  method  employs  techniques  in  non-linear,  least  squares  approx. mat  ions  and  is  applicable  to 
both  high-  and  low-order  transfer  functions  (•«).  On  the  final  Interpol 1 1 Ion  in  the  program,  the  least 
squares  error,  LSE,  of  the  best-fitting  transfer  function  is  obtained  as 

L13E  *=  YTY,  (b) 

where  Y  is  a  residuals  matrix  (Y^i  transpose  of  Y)  which  is  made  up  of  error  terms  that  express  the 
differences  between  the  experimental  data  and  the  model  which  is  to  be  fitted.  The  mean  square  error, 

m:;e  -  LSE/(2L  -  S),  (?) 

compares  the  goodness-of-fit  of  the  derived  transfer  function  to  that  of  the  experimental  data  (2L  f 
number  of  real  and  imaginary  components  of  the  data  points,  and  i  number  of  parameters  in  the  transfer 
function). 

A  more  detailed  description  of  the  methodology  may  be  found  in  a  special  monograph  (7)  anc  elsewhere 

(6,  21). 


RESUIX-  AND  DISCUSSION 

The  simplest  transfer  function  that  fitted  the  data  fur  all  units  is  of  the  form. 
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where  s  i it  tract  ion-.l  —  rd«'r  differential  i>per»t"r  with  0  <  k  <  1,  and  0  is  a  gain  constant  with 
units  in  impulut-::«s”V<i*  grct-s  •:»*■•*  (-’1).  The  Bode  plntt;  for  four  of  the  units  are  shown  in  Figs.  2-5, 
together  with  their  best -fitting  »'(:;),  and  the  best-fitting  tors  ion -pendulum  model 


d-Ul  * 


(9) 


(In  rigs.  2-5,  the  effects  of  the  x.?  -  i’.O  ms  (vide  infra)  contributes  no  more  than  to  the  gain  spectra 
between  f  =  0.01  and  10  Hz;  consequently ,  the  single-pole  transfer  function,  G"(s),  was  an  adequate 
representation  of  the  torsion-pendulum  model.)  As  is  evident  ip  the  plots  in  Figs.  2-5,  G'(s)  is  a 
much  be*  ter  *'it  to  the  lata  than  is  G"(n)  (of.  the  MSF.  values  for  the  two  models);  and,  in  particular, 
t he  fit  appears  to  improve  with  increasing  k. 

Wnat,  then,  is  the  significance  of  k,  or  better  still,  s  ?  Elsewhere,  it  is  shown  that 


M  D4  s  U.^s+l) 
i=l  ( i ,s+l  ) 


(10) 


where  K  and  I>j  are  constants,  M  *  •  (in  theory,  but  finite  when  fitting  Eq.  (10)  via  a  digital 
computer),  and  and  ( C0  i  0)  are  time  constants  (21).  Interestingly,  when  M  =  1,  Eq.  (10) 

becomes 


(11) 


which  has  previously  been  defined  as  the  transfer  function  of  the  adaptation  operator  (27,  bo).  Thus, 
s^  appears  to  be  a  form  of  adaptation.  Work  by  Thor son  and  Biederman -Thor son  (39)  suggests  that 
s11  represents  a  distributed  relaxation  process  which  is  inherent  in  the  sensory-adaptation  mechanics 
of  Li mu 1  us  photoreceptors,  vertebrate  retinal  receptors,  chemoreceptors ,  and  other  mechanoreceptors . 
Investigations  by  Taglietti,  Hossi  and  Casella  (38)  further  suggest  that  s*  likely  represents  a 
relaxation  phenomenon  consisting  of  a  time-varying  intracellular  electrogenic  process,  the  com¬ 
ponents  of  which  are  summed  with  the  generator  potential  In  the  receptor  hair  cell. 


The  coefficient,  of  variation,  CV,  was  determined  as  the  ratio  of  the  standard  deviation  of  intervals 
to  the  mean  interval,  as  obtained  from  int erspike-Interval  distributions  of  spontaneous  single  unit 
activity.  When  a  regression  of  CV  on  k  was  made  for  2 8  units,  a  statistically-significant  product- 
moment  correlation  (r  *  0.3ftli,  F»  <  0.05)  was  obtained  (21).  Thus,  the  larger  the  CV  is,  the  larger  the 
value  of  k  and,  consequently,  the  amount  of  adaptation.  Other  work  by  Goldberg  and  Fernandez  (15).  in 
squirrel  monkeys,  shows  that  the  CV  is  statistically  correlated  with  semicircular-canal  afferent 
fiber  conduction  rates.  The  thicker  fibers  have  faster  conduction  rates  and  larger  CVs.  Together, 
these  findings  suggest  that  sensor >  adaptation  phenomena  are  directly  dependent  on  the  innervation 
pattern  of  the  afferent  fibers. 

The  transfer  function,  G*(s),  differs  also  from  that  of  G”(s)  in  that  is  not  single-valued  as 
it  is  in  the  torsion-pendulum  model;  rather,  it  is  unit  dependent,  taking  on  values  from  x^  *  U.U5  to  22.17 
a  (mean  ±  SEM  -  10. 2U  *  1.20  s)  (21).  (In  fitting  Eq.  (ft),  the  coupling  between  k  and  x^  would  account 
for  some  of  the  five-fold  range  of  values  that  were  determined  for  tL.  However,  there  is  sufficient 
indication  from  other  studies  (32)  that  the  response  dynamics  of  small  groups  of  contiguous  hair  cells 
are  quite  different  from  those  of  other  groups.)  Kealizing  that  the  nair-cell  tufts  are  stiff  (lb), 
that  their  lengths  vary  according  to  their  position  on  the  sensory  epithelium  (crista)  (22),  that  the 
number  and  thickness  of  the  stereocilia  can  be  variable  (22),  and  that  the  mechanical  properties  of  the 
cupula  are  not  necessarily  uniform  across  the  crista  (29),  then  it  is  plausible  that  x^  could  have 
a  regional  distribution. 


The  form 

.  _  Cb* 

‘  (TIB*l)(T..Sn)  •  (12) 

or  sometimes, 

0(s)(t  s*1)  (13) 

m 

was  fitted  to  published  afferent-response  data  in  the  squirrel  monkey  (Salmirl  sciureua)  (12),  the 
frog  ( Ran a  esculenta  and  R.  temporarla)( 3) .  the  gerbil  (Merlones  ungulculatusT  ( 35 ) ,  and  the  gult&rflsh 
(Rhlnobatos  productus)  (33).  The  parameter  xm  =  high  freauency  time  constant  which  results  from  both 
the  displacement  and  the  rate  of  displacement  of  the  cupula.  Such  a  term  has  been  obtained  from 
analysis  of  vestibular -driven  eye  movements  in  man  (Benson  and  Sternfeld,  cited  in  ref.  l),  and 
primary  afferent  canal  responses  in  the  squirrel  monkey  (12)  and  the  elasmobranch  fish  (2l*).  Table  1 
lists  k,  t^,  t t. ,  and  for  these  four  species  and  for  unit  responses  to  white-noise  stimuli  obtained 
from  pigeons  that  were  primarily  enc^phale  isoi£  preparations.  For  all  species  listed  in  Table  1,  the 
MCE  using  Eqs.  (12)  or  (13)  was  comparable  to  or  smaller  than  that  obtained  with  other  models.  Further 
details  are  given  elsewhere  (8). 
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:  Parameters  >.  f  a  genera. i.* 
dynamics  in  five  se 

?  .not.  i  .ti  (Eqs.  ( 

i.')  or  (13))  describing  semi ci rcu J ar- 

7 

- T 

t , 

, 

1 requenoy  range , 

Species 

4  i  ms  I 

m 

(s> 

f  !\\7  ) 

Squirrel  mohrtey 

(a)  "regular”  uni*. 

i 

- 

0.0125-8 

(b)  "irregular"  u. . * 

1 

0.03 

0.0125-8 

Gerbi 1 

(a)  cv  <  ■:  1 

2 

- 

0.01-5 

(b)  fV  >  0.1 

.*  ; 

- 

0.01-5 

Frog 

- 

0.0125-0.5 

Guitar fish 

" 

0.U2-L 

Pigeon 

.  I.  2 

- 

0.5-16 

Table  l,  the  frequency  » 

.  .  t  restricted  in 

the  frog  and 

guitar  fish  to  utilize  a 

adequate  value  of  r<;.  In  genei  n 
from  the  biophysical  proper* ie- 
of  the  membranous  vestibular  a;  i  *> 
Money  and  colleagues  ( U )  f  r  *• 
(f  =  of  80  Hz  in  the  rv- 

from  canal  afferents  in  the  g  : 
given  evidence  that  5.70  ^  t  • 
fitting  to  afferent  data  fn  t  h«* 


(=  0/H)  in  the  fits  have  been  predetermined  indirectly 
ymph  and  the  dimensions  of  the  pertinent  anatomical  features 
.•u  ur  trj  =  2  ms,  which  was  determined  biophysical  1  y  by 
.vr.it;  that  there  should  be  an  upper  break  frequency 
«  P  i.  >e  dynamics  in  semicircular-oanal  afferents.  Recording 
•  ..•is  frequency  up  to  70  Hz,  Hartmann  and  Klinke  ( >  have 
1 1,  the  same  range  as  those  determined  empirically  for 
.  ,  .'1 ,  45). 
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